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ABSTRACT 

We propose a model to account for the bolometric light curve, quasi-continuum 
and the Ca ii emission features of the peculiar type la supernova (SN) 2002ic, which 
exploded in a dense circumstellar envelope. The model suggests that the SN la had 
the maximum possible kinetic energy and that the ejecta expand in an approximately 
spherically symmetric (possibly clumpy) circumstellar environment. The Ca ii and 
quasi-continuum are emitted by shocked SN ejecta that underwent deformation and 
fragmentation in the intershock layer. Modeling of the Ca ii triplet implies that the 
contribution of the O i 8446 A line is about 25% of the 8500 A feature on day 234, 
which permits us to recover the flux in the Ca ii 8579 A triplet from the flux of 8500 
A blend reported by Deng et al. (2004). We use the Ca ii doublet and triplet fluxes on 
day 234 to derive the electron temperature {~ 4400 K) in the Ca il line-emitting zone 
and the ratio of the total area of dense fragments to the area of the shell, S/ Sq ~ 10^. 
We argue that Ca II bands and quasi-continuum originate from different zones of the 
shocked ejecta that reflect the abundance stratiflcation of the supernova. 
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1 INTRODUCTION 

Recently, strong evidence has been presented that a sub- 
class of type Tin supernovae (SN) ('n' stands for narrow Hq) 
includes rare SN la events exploding in a dense circumstel- 
lar medium (CSM) (Hamuy et al. 2003). To date, the new 
variety comprises three events: SN 2002ic, SN 1997cy and 
SN 1999E (cf. Deng et al. 2004). Statistical arguments indi- 
cate that the new family constitutes less than one percent 
of all SN la (Chugai & Yungelson 2004). All these super- 
novae remain very luminous after the light curve maximum, 
which is related to the CS interaction. A light curve model of 
SN 1997cy suggests that the CS envelope within 3 x 10^* cm 
contains several solar masses (Chugai & Yungelson 2004). 

Two options for the origin of the CSM have been pro- 
posed (Hamuy et al. 2003): mass-loss by a supergiant in a 
binary scenario of SN la (Whelan & Iben 1973) or the mass- 
loss by a single supergiant in the SN1.5 scenario (Iben & 
Renzini 1983). The similarity of SN 2002ic events seems to 
favour the SN 1.5 scenario, with the progenitor mass around 
8 Mq (Chugai & Yungelson 2004). The possible scenario of 
a C/O white dwarf (WD) merging with the CO-core of a 
supergiant companion (Livio & Riess 2003) faces a serious 
problem in explaining the short time lag between the ejec- 



tion of the common envelope and the explosion (Chugai & 
Yungelson 2004). 

The study of the CS interaction and modeling of the 
spectra should elucidate the late evolution and the origin of 
these interesting events. An important recent conclusion is 
that SN 2002ic interacts with an equatorially concentrated 
CSM (Deng et al. 2004). This statement is based upon 
the interpretation that line profiles in SN 2002ic are too 
broad compared with the predictions of a spherical model. 
Also, spectropolarimetric observations of SN 2002ic reveal 
the presence of polarization which is likely to be related to 
the CSM (Wang et al. 2004). The claimed asymmetry of the 
CSM may have an important implication for the problem 
of the evolutionary scenario, because the single and binary 
scenarios are likely to possess difi'erent degrees of symmetry. 

The aspherical model predicts variations of the emis- 
sion line profiles of SN 2002ic-like supernovae viewed from 
different angles. However, all three known objects (see Deng 
et al. 2004, their Fig. 1) show similar shapes and widths 
of spectral features. This fact implies that, despite the de- 
tected polarization of the emission from SN 2002ic, the in- 
teraction possesses a high degree of sphericity. In view of 
this, we investigate a spherical model here. We attempt 
to explain the light curve, the Ca ll line profiles and the 
quasi-continuum of SN 2002ic-like events on the basis of a 
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circumstellar interaction model. We rely on the spherical 
approximation, although we allow for local deviations from 
sphericity (e.g. clumpy structure). We start with the simula- 
tion of the CS interaction dynamics and the bolometric light 
curve of SN 2002ic using the thin shell approximation (Sec- 
tion 2). Although the model is the same £is before (Chugai 
& Yungelson 2004), we include more realistic ejecta abun- 
dances, which affects the light curve. We then analyse the 
line formation in the inhomogeneous intershock layer and 
propose a simple treatment of this problem (Section 3). Us- 
ing this approximation we simulate the Ca ll line profiles 
and the quasi-continuum of SN 2002ic with the kinematic 
parameters of the interaction model. We demonstrate that 
the spherical model successfully accounts for the basic ob- 
servational properties of SN 2002ic and discuss implications. 



2 CS INTERACTION AND THE LIGHT 

CURVE 

The interaction of a SN la with a dense spherical CS enve- 
lope and the light curve formation is calculated using the 
model applied earlier for the analysis of the CS interaction 
of SN 2002ic (Chugai & Yungelson 2004). We recapitulate 
the basic assumptions of the model and emphasize some 
modifications. 



2.1 Model 

The density distribution in a freely expanding SN la ejecta 
{v oc r) is approximated by the exponential law p oc 
exp(— w/wo) with vq determined by the the mass M = 
1.4 Mq and kinetic energy E. An exponential law better 
represents the density distribution in SN la ejecta than a 
power law (Dwarkadas & Chevalier 1998). In an attempt 
to provide the maximum expansion velocity of the interac- 
tion shell we consider a SN la model with the maximum ki- 
netic energy, corresponding to the detonation model DETl 
of Khokhlov et al. (1993): E = 1.75 x lO" erg with the 
^®Ni mass of Mni = 0.92 Mq. We do not claim that the 
DETl detonation regime actually occurs in SN 2002ic. This 
model strongly underproduces intermediate mass elements 
(e.g. Ca) (Khokhlov et al. 1993), which is at odds with the 
strong Ca n lines observed in SN 2002ic (see below Sec. 
3.4). Parameters of the DETl model are adopted because 
this case illustrates what happens if all of a Chandrasekhar 
mass C/0 WD is burned. The important point here is that 
the complete incineration of C and O is not required: if the 
burning in the outer layers is terminated at the intermediate 
elements, the released rmclcar energy is practically the same 
as in the case of complete incineration. 

The CS density distribution is described by power laws 
(p oc r'°) with different power law indices for each of three 
suggested radial zones: ki in the range r < ri, k-z for 
ri < r < r2, and fes for r > r2. The variation of the power 
law index is suggested by two observations: (i) the luminos- 
ity decline after about day 400 (Deng et al. 2004) suggests 
a CS density drop beyond some radius r > r2; (ii) the ris- 
ing early visual luminosity of SN 2002ic (Wood-Vasey et al. 
2002; Hamuy et al. 2003) suggests a flat (or rising outwards) 
CS density distribution in the inner part (r < ri ) of the CS 
envelope (Chugai & Yungelson 2004). 
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Figure 1. Bolometric light curve (upper panels), radius of the 
thin shell (middle), and expansion velocity (bottom) for models 
with CS density p oc r^^ (a) and p oc r^^'^ (b). The bolometric 
light curve with CS interaction {solid line) and without CS inter- 
action {dotted) arc ovcrijlottcd on tiic data of Deng ct al. (2004). 
The lower panel shows the velocity of the thin shell {solid) and 
the SN boundary velocity {dashed). 



The hydrodynamic interaction is calculated in the thin 
shell approximation, in which the layer between the forward 
and reverse shocks is replaced by a geometrically thin shell 
(Chevalier 1982b). The kinetic luminosities of the forward 
and reverse shock waves are converted into X-ray luminosi- 
ties using transformation coefiicients dependent on cooling 
rates of the postshock gas and the expansion time (Chugai 
1992; Chevalier & Pransson 1994). For the CS density range 
of interest, the reverse shock is always radiative, which leads 
to the formation of a cool dense shell (CDS) composed of the 
shocked SN ejecta. The forward shock can be also radiative 
and contribute to the CDS formation. The X-ray radiation 
of both shocks absorbed by the cool material (SN ejecta, 
CDS, and CSM) is identified with the optical bolometric 
luminosity powered by the CS interaction. The additional 
luminosity powered by the radioactive decay *®Ni - ^®Co - 
^^Fe is calculated using an analytical theory for the SN light 
curve (Arnctt 1980). The resulting light curve is the sum of 
the interaction and radioactive luminosities. 

An important modification in the present version is the 
inclusion of a 'realistic' chemical composition of the SN la 
ejecta instead of the previously assumed solar composition 
(Chugai & Yungelson 2004). The SN la composition is now 
approximated by a mixture of Fc and Si, where the amount 
of Fe is equal to the ejected ^^Ni mass. The new composition 
results in more efficient cooling in the reverse shock and 
more efficient absorption of the X-rays in the SN ejecta. The 
adopted absorption coefficient for X-rays in the SN material 
is kx fti SOOS^^y cm^ g~^, i.e., 50 times larger than that 
with a solar composition. Owing to the efficient absorption 
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Figure 2. Circumstellar density in tlie model witli p oc [thin 
line) and p oc * {thick line). In both models, the inner zone 
has a plateau density distribution to provide the rise of the light 
curve to the maximum at about day 20. 

of X-rays in the SN la ejecta with the 'reaUstic' composition, 
the rate of emission of optical photons per unit of kinetic 
luminosity is higher in the present model compared with 
the earlier version. 



2.2 Light curve 

Our modeling of the bolometric optical light curve is com- 
pared with the SN 2002ic data from Deng et al. (2004), mod- 
ified here for a Hubble constant of Ho = 70 km s~^ Mpc~^. 
Two versions of the light curve computations along with 
the corresponding shell radius and velocity are presented in 
Fig. 1. Fig. la shows the model with k2 = —2, whereas the 
second model (panel b) corresponds to k2 — —1.4. Both CS 
density distributions are shown in Fig. 2. The model with 
k2 = —1.4 provides the largest expansion velocities of the 
thin shell (vs) and of the SN at the boundary (vsn). From 
the point of view of the interpretation of the observed broad 
line profiles, this model is preferred. Henceforth, this case is 
referred to as the standard model. 

The flatter CS density distribution of the model com- 
pared with a steady-state wind implies that after the main 
mass loss episode, the mass loss rate gradually decreased 
while the presupernova evolved towards the explosion. The 
total mass of the CS envelope in the standard model is 
1.6 Mq. This estimate should be considered a lower limit, 
because we neglect the possible wind at r r2 that may 
comprise a significant mass, despite its lower density. If the 
wind velocity out to r2 ~ 2.6 x 10^^ cm is > 10 km s~^, then 
the age of the strong mass loss is < 820 yr and the mass loss 
rate is > 2 x 10"^ M© yr"^ 

We recalculated the light curve of SN 1997cy for the 
detonation model with a 'realistic' composition and found 
that the mass of the CS envelope in this case is « 2.5 Mq, a 
factor two lower than our previous estimate for SN 1997cy 
based upon a model with lower explosion energy and normal 
composition of the SN ejecta (Chugai & Yungelson 2004). 
Most of the difference is related to the efficient conversion of 
the kinetic luminosity into the bolometric luminosity owing 
to the efficient absorption of X-rays by the ejecta. 



The model of the optical light curve predicts the lumi- 
nosity and flux of the escaping X-rays, which may be used 
as an additional test of the model for SN 2002ic-like events. 
In Fig. 3, we plot the total luminosity of the escaping X- 
rays from the forward and reverse shocks, together with the 
temperatures of both shocks for the standard model. Due 
to the high CS density and rather flat CS density distribu- 
tion, the velocity jump at the inner shock is relatively high 
(Fig.l). This, combined with the higher molecular weight 
of the SN la matter, results in a high temperature of the 
reverse shock that is comparable to the temperature of the 
forward shock after day 250. In a SN II interacting with a 
wind, the reverse shock is usually markedly slower and cooler 
than the forward shock (Chevalier 1982b). Remarkably, the 
luminosity of the forward shock is large and provides a sub- 
stantial contribution to the optical SN luminosity. This is 
due to the fact that the forward shock is radiative: it enters 
the fully radiative regime on day ~ 150. The reverse shock 
is radiative through the time period we have modeled. The 
inflection in the forward shock luminosity on day ~ 400 re- 
flects the corresponding inflection in the adopted CS density 
distribution at r « 2.6 x 10^® cm (Fig. 2). The forward shock 
becomes nonradiative soon after that point, at f > 440 d. 

The photon spectrum of the emergent X-ray flux 
at different epochs between day 30 and day 600 is 
shown in Fig. 4 together with the detection limit of the 
JEM-X monitor of the INTEGRAL satellite assuming a 
3a detection level with the integration time of 10" s 
(http://astro.estec.esa.nl/Integral/). The plot demonstrates 
that the detection of X-rays from SN 2002ic is not expected 
in our model for a reasonable integration time. On the 
other hand, a SN 2002ic-like event exploding at a distance 
D < 100 Mpc would be detectable at the favorable epoch 
of 200-400 days by the JEM-X monitor with an integration 
time < 10® s. Chandra and XMM may be even more efficient 
in detecting the X-rays from distant SN 2002ic-like events. 
Between days 200-500 the predicted flux of SN 2002ic in the 
0.5 — 10 keV band is ~ 10^^* erg s^^ cm^'^. Given the Chan- 
dra detection limit ~ 10~^^ erg s~^ cm~^ in the 0.5—10 keV 
band, we conclude that X-rays from SN 2002ic could have 
been detected at an age 200 — 500 days. At later epochs, the 
flux quickly decreases (Fig. 4) and on day 600 it becomes 
far below the Chandra sensitivity. 

The above model for the escaping X-ray flux assumes a 
smooth GSM. dumpiness may modify these results in two 
ways: a different temperature in the forward shock, and a 
different escape probability for photons through the clumpy 
GSM. 



3 FORMATION OF THE OPTICAL 

SPECTRUM 

The late time Subaru spectrum on day 234 (assuming the ex- 
plosion date JD 2452585) consists of broad emission bumps 
(we dub them quasi-continuum) with superimposed emission 
lines of Ha, H/3, Ca ii 3900 A, and Ca ii 8600 A blended with 
O I 8446 A (Deng ct al. 2004). Most of the bolometric lu- 
minosity in our dynamical model comes from the shocked 
ejecta. We, therefore, believe that, apart from Ha and H/3 
and other narrow lines, most of the optical spectrum on day 
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Figure 3. Luminosity of escaping X-rays (top) and tempera- 
ture (bottom) of the reverse shock (solid line) and forward shock 
(dashed) 



234 (i.e. quasi-coritirmum and Ca ll linos) originates from 
shocked SN ejecta. This conjecture will be checked below. 

We only briefly address the important issue of Ha line 
formation. Although the CS origin of the narrow Ha emis- 
sion is certain, the nature of the broad Ha is unclear. Mul- 
tiple Thomson scattering as a mechanism for the broad Ha, 
applied earlier to SN 1998S (Chugai 2001), and discussed by 
Wang et al. (2004) in connection with SN 2002ic, requires 
a large optical depth tt ~ 3. This is in conflict with our 
dynamical model that predicts a Thomson optical depth of 
the undisturbed CS ejecta of tt ~ 0.1 on day « 240 . In 
our opinion, the broad Ha comes from the shocked cool 
dense hydrogen. There are three possible origins for that 
line-emitting gas: (i) the outer layers of the SN ejecta may 
contain the remains of the hydrogen envelope of the prcsu- 
pernova (in the SN 1.5 scenario), so the cool shocked ejecta 
enriched by hydrogen may be partially responsible for the 
broad Ha component; (ii) the forward shock becomes rar 
diative in the case of smooth CS hydrogen after about day 
150, so the shocked CS hydrogen cools at that epoch and 
becomes able to emit broad Ha; (iii) in the case of a clumpy 
CS envelope radiative shocks in CS clouds and the shear flow 
around the shocked clouds may produce the broad velocity 
spectrum of the line-emitting gas that can be responsible 
for the intermediate and broad Ha line. In reality, all three 
sources of the cool shocked hydrogen may contribute to the 
broad Ha component. However, only in the third option can 
the observed intermediate Ha component with a width of 
~ 2 X 10^ km s"^ (Hamuy et al. 2003) be produced. 

Henceforth, we concentrate on the study of the Ca ii 
emission lines and the quasi-continuum. 



3.1 Line formation model 

Momentum and mass conservation imply a very small thick- 
ness of the CDS, of the order of 5 ~ i?(vsh/cs)^, where 
Vsh = Vsn — Vs is the velocity of the reverse shock and Cs 



is the ion sound speed in the cool gas. With Vsh ~ 3000 km 
s~^ at an age of 230-250 d (see Fig. 1), adopting the tem- 
perature of the cool gas of lO** K and the ion average atomic 
weight ^ 40 (i.e., Cs « 2 km s"^), 5/R ~ 10"®. For such a 
thin shell the velocity dispersion along the line of sight is 
negligible, so the line radiation transfer in this shell may be 
treated in the static approximation (Chevalier and Pransson 
1994). The expansion effect appears only in the transformer 
tion of the shell emission to the observer frame. However, 
the picture of a spherical thin CDS is an idealization; in 
reality, as we will show below, the shell should be strongly 
distorted. 



3.1.1 Line-emitting shell 

A straightforward argument against the existence of a ge- 
ometrically thin spherical shell is based upon the observed 
profiles of the Ca ii features. It is well known that in the 
optically thin case the line profile from a geometrically thin 
expanding shell should be boxy. This might be a sensible 
approximation for the observed Ca ll features given some 
unavoidable perturbation and smearing of the boxy profile. 
The problem, however, is that the Ca II lines are expected 
to be optically thick, not thin. 

To show that the dense shell is opaque in the Ca ii 
linos, wo consider the conditions on day 240. The swept up 
SN mass in the reverse shock in the standard model is ~ 
0.6 Mq and the shell radius is _R ~ 1.8 x 10^^ cm. Assuming 
that the swept up mass resides in the CDS and that the Ca 
abundance among Fc-jjeak elements is solar, i.e., « 0.03 by 
mass (a sensible approximation for a rough estimate) one 
gets the column density of Ca ii in the CDS of Ai'(Call) ~ 
2 X 10^° 02 (where cj)2 is the ionization fraction of Ca ii). The 
optical depth in the Ca ii doublet lines is then 



^^iV(CaII)~10% 
mc ut 



(1) 



where /12 = 1 is the oscillator strength for the Ca 11 doublet 
and Mt ~ 2 km s~^ is the thermal velocity assuming T = 
10* K. On day 234, the number density in the CDS under 
pressure equilibrium is lO^"* cm^'*. Assuming LTE in the 
cool dense gas we obtain 10"^ < 02 < 1 for 4000 < T < lO"* 
K. The equation (1) shows that the optical depth of the CDS 
in the Ca 11 doublet is very large for the estimated Ca 11 
ionization fraction. It is easy to check that, given the low 
excitation potential of the Ca 11 term ^D, the Ca 11 triplet 
should be optically thick too. 

Optically thick Ca 11 lines from a geometrically thin 
shell are expected to be double-horned (M-shaped) F cx 
A — A()| with peaks at the shell velocity (±Ws) (Gerasimovic 
1933). This type of profile is apparently inconsistent with 
the smooth, round-topped (flat-topped) Ca 11 doublet and 
triplet proflles observed in SN 2002ic (see Deng et al. 2004). 
To resolve this controversy, we follow Cid Fernandes & Ter- 
levich (1994) in suggesting that the round-topped Ha line 
proflle observed in SNe II interacting with a dense CS envi- 
ronment arc produced by a thin shell corrugated due to the 
Ray Icigli- Taylor (RT) instability. 

It is well known that the decelerated thin shell in the 
intershock layer is subject to the Rayleigh- Taylor (RT) in- 
stability (Chevalier 1982a; Nadyozhin 1985). RT spikes grow 
and degrade via shearing and advective flows, thus produc- 
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Figure 4. X-ray photon flux from SN 2002ic predicted by the 
standard model at different epochs between days 30 and 600 as- 
suming Ho = 70 km s^^ Mpc^^ {D = 285 Mpc). The 3a sensitiv- 
ity curve of JEM-X monitor of INTEGRAL (dashed line) assumes 
a 10® s integration time. 



ing a two-fluid mixture consisting of dense cool gas and hot 
postshock gas behind the forward shock wave. As a result, a 
mixing layer with an extent AR/R ^0.1 forms between the 
thin shell and the forward shock as demonstrated by 2D and 
3D simulations (Chevalier & Blondin 1995; Blondin & El- 
lison 2001). The clumpy structure of the CSM may cause 
additional perturbations of the mixing layer and lead to 
the corrugation of the forward and reverse shocks (Blondin 
2001, his Fig. 6). If both reverse and forward shocks are 
radiative, the shell instability becomes even more vigorous 
due to the nonlinear thin shell instability (Blondin & Marks 
1996; Hueckstaedt 2003). The above arguments beised upon 
the profile shape and the interaction hydrodynamics suggest 
that a realistic model of the intershock line-emitting layer 
should include the inhomogeneous distribution of the dense 
shocked gas. 

To treat the line formation in the inhomogeneous inter- 
shock layer, we reduce it to a spherical 'shell' of a thickness 
AR = 2^R = O.IR in the range of radii \r - R\ < ^R 
filled with similar, randomly oriented, dense planar frag- 
ments (Fig. 5). Each fragment is a thin sheet with an area 
a and total surface area 2a. As a result of the deformation 
and stretching of dense fragments, their cumulative area S 
should, generally, exceed the area of the unperturbed spher- 
ical dense shell So ~ ^tvR^ by a factor S/So > 1. As we will 
see below, this parameter (called the 'area ratio') turns out 
to be crucial for the interpretation of the high luminosity of 
the Ca II features in SN 2002ic. Moreover, the high lumi- 
nosity of the Ca ll features requires S/So ^ 1. This is the 
reason why we omit in our simple model (Fig. 5) the thin 
spherical shell that is produced and destroyed continuously 
in the reverse shock. 



3.1.2 Line radiation transfer in the shell 

We consider two plausible velocity distributions of filaments 
in the shell R{1 — f) < r < R{1 + ^): a constant veloc- 
ity and a linear law v = Vs{r/R). Line photons emitted by 
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Figure 5. A schematic view of SN 2002ic that shows the frag- 
mented cool dense shell resijoiisiblc for the bulk of the optical 
radiation. The fragmented shell is bounded by the forward shock 
propagating in the CSM and by the reverse shock propagating in 
the unshocked SN ejecta. 



randomly distributed dense fragments may escape or be ab- 
sorbed by another fragment of the shell. To describe the line 
radiative transfer inside the inhomogeneous shell we apply 
the local escape probability language adapted for the 'gas' 
ofj opaque macroscopic chunks of dense material. Definitions 
of the local optical depth and emissivity must be modified 
accordingly. 

The total number of fragments in the shell is A'^ = 
AivR? {S / So) / a and the fragment number density is rif — 
N/V = {S/So)/{aAR), where V = A-kR^AR. The average 
cross-section of a fragment is (cr) with the ratio {a) /a = 0.5 
for plane sheets and ~ 0.5 for deformed sheets. For instance, 
for a cylindrical surface with a radius a and a height h <^ a 
this ratio is ~ 0.5 with the accuracy of 0{h/a). The radi- 
ation escapes from an optically thick fragment in the fre- 
quency band 2xAi/o assuming a gaussian absorption coeffi- 
cient fc,y oc cxp [— (A/z/Ai/d)^], where Auu ~ 'f-t/A and x is 
defined by the condition t{x) = 1, i.e., x = (In r)" '\ A frag- 
ment absorbs the incident radiation in the same frequency 
band, 2xAi/d. 

The analog of the Sobolev optical depth in our 
case is the local occultation optical depth, Q — 
2xut{ij)n(\dv/ ds\^^ . To take into account the finite opti- 
cal depth of a fragment this expression should be multiplied 
by the absorption probability of a fragment (1 — e~'^). For 
a shell with constant velocity, the velocity gradient along 
some direction is dv/ds = (vs/R) sin^ 6, whore 9 is the an- 
gle between the photon wave vector and the radius. The 
expression for Q then reads 



Q 



(I) 



\ARj Us 
= Qo(l-e 



(1-e- 
^)/sin'( 



(2) 



where (o")/a" = 0.5 is used. For the linear velocity case, 
Q = Qo{l — e~^). To estimate Qo let us adopt Vs = 6000 
km s~^, Wt = 2 km and a; « 3. Equation (2) gives 
Qo ~ 10~^(S/So). Therefore, the shell becomes optically 
thick (Qo > 1) for an area ratio S/Sa > 10^ provided r 3> 1. 

The emissivity of the fragment ensemble is j = 
(l/27r)crnfF, where F is the frequency integrated flux. In 
the isothermal approximation and assuming Boltzman pop- 
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8500 A is dominated by O I 8446 A as Deng at al. (2004) 
suggest. Below we will see that the contribution of the O I 
8446 A line is small so the high luminosity of the Ca ll triplet 
is primarily due to the large area ratio. 

To show the effect of Qq wc consider the transparent 
case Qo = and modestly opaque shell Qo = 1 for t 3> 1. 
The calculated profiles for a constant velocity shell with 
Vb = 6000 km s~^ and ^ = 0.05 are shown in Fig. 6 (upper 
panel). The finite optical depth smooths the boxy profile as 
expected. In the limit of Qo 1, the profile is parabolic 
with a weak asymmetry due to relativistic effects. The rel- 
ativistic effects arising from the Lorentz transformation of 
the frequency, wave vector and intensity arc usually ignored 
in SN spectrum calculations. However as was demonstrated 
by Jeffery (1993) these effects become noticeable even for 
expansion velocities of ~ lO** km s~^ and generally result in 
the skewing of an emission component towards blue. 



Figure 6. The line profile for two different models of the cold, 
dense shell. Upper panels show profiles for a transparent {Qo = 0) 
and a modestly opaque shell (Qo = 1). Lower panels show profiles 
for the same models, but with the addition of non-local scattering 
in the SN ejecta for a cut-oflF radius rc = 0.95 {thick line) and 
rc =0.9 {thin line). 

ulations, the fiux is F = 2itB,^{T)xAi/o{1 — e""^), so the 
emissivity is 

j = anfB^{T)xAuD{l-e~''). (3) 

To compute the luminosity, one needs to take into ac- 
count self- absorption via the angle-averaged escape proba- 
bility /3 = ([1 — exp(— Q)]/Q). Assuming pure absorption 
one gets the luminosity L = AirjVp, where V is the shell 
volume AitB^AR. With the line emissivity j from equation 
(3), the luminosity of Ca ii triplet is 

L = lQ-K^R^{S/So)B^{T){\ - e-^)PxAuu . (4) 

The luminosity increases linearly with the area ratio for 
Qo 1 (i.e., /3 « 1) and saturates for large S/So, when 
Qo becomes also large. Actually, for large Qo one gets 
(3 ~ 1/Qo oc So/S so L becomes independent of S/So- 
Remarkably, the line luminosity according to equation (4) 
exceeds the upper limit on the luminosity of a static shell 
by a factor {S/So)f3 that, for large area ratio, may be sub- 
stantially greater than unity. This apparent paradox is re- 
lated to the velocity gradient that permits the volume emis- 
sion regime for resonance photons as opposed to the surface 
emission regime in the static case. A finite albedo of a 
fragment for incident photons requires a modification of the 
escape probability: /3 in equation (4) should be replaced by 
p = /3/[/3+{l—f3){l—u))]. However, we compute the line pro- 
files on the assumption of zero albedo, which, as we checked, 
does not affect the profile significantly. 

In the limit of a small area ratio (5/ So ~ 1) the lumi- 
nosity of the Ca ll triplet on day 234 for the probable tem- 
perature range T < lO"* K for 5 = 1 is L(8579A) < 5 x 10''° 
erg s~^. On the other hand, using the fiux in the 8500 A 
band, / ^ 4.7 x 10"^" erg s"^ cm"^ on day 234 (Deng et al. 
2004), we find the luminosity L « 3.8 x lO"^ erg s"^ which 
significantly exceeds the upper limit found for S/So = 1. 
This means that either S/So ^ 1 or the emission feature at 



3.1.3 Non-local scattering in SN ejecta 

Although the contribution of the net emission from the un- 
shocked SN ejecta is small for SN2002ic given the domina- 
tion of the shock luminosity, the SN ejecta can affect the 
line profile via the scattering of line photon omitted by the 
CDS. If the SN ejecta are optically thick in the same line, 
then a line photon emitted by the shell inwards can scatter 
on the unshocked Ca ii ions at the resonance point with the 
distance from the emission point 

s = RM (l-i^) , (5) 

where fj, = cos ^ < 0. The resonance points reside on a 
sphere of radius (1 — Vs/vsn)R/2 that contains the emitting 
point (s = for /i = 0). 

To compute the effect of the non-local scattering of line 
photons in SN ejecta we adopt the boundary velocity of the 
SN ejecta Vsn = 9000 km s~^; this value and the shell ve- 
locity (6000 km s~^) are close to velocities in the standard 
model around day 230 (Fig. 1). Let rc bo a cut-off radius (in 
units of Rs) defined in such a way that the Sobolev optical 
depth of the SN ejecta Tsn = in the outer zone, r > Vc, 
and Tsn > in the inner zone, r < rc. In both models the ef- 
fect of non-local resonance scattering is shown for two cases: 
Tc ~ 0.95 and rc = 0.9, assuming Tsn = 5 in the unshocked 
SN (Fig. 6, lower panel). Profiles demonstrate significant 
skewing towards the blue and the emergence of a symmetric 
broad base produced by scattered photons in the fast SN 
ejecta with velocity Vsn > Vs. The shape of the line profile 
depends on rc: the skewing is larger for larger rc. As we see 
below, the non-local scattering can account for the asymme- 
try of the Ca n doublet observed in SN 2002ic. 

3.2 Models of the Ca ii features 

We start with the modeling of the Ca ii doublet. A Monte 
Carlo technique is used for the radiation transfer computa- 
tions. The parametrized emission rate of the CDS fragments 
in linos of the Ca ll doublet is set assuming oquipartition of 
the upper levels. The non-local interaction in the shell is 
treated as a pure absorption process, whereas the resonance 
interaction of line photons with the unshocked SN ejecta is 
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Table 1. Model paramet 
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0.95 
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DCl 


1 


5936 


9092 


0.95 
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DLO 





5936 


9092 


0.95 


3 




DLl 
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5936 


9092 


0.95 


3 




TLO 





6053 


9230 


0.95 


0.005 


0.25 


TLO 


1 


6053 


9230 


0.95 


0.005 


0.25 



considered as pure scattering. The quasi-continuum radia- 
tion, presumably emitted by the shell cam be scattered by 
Ca II lines of SN ejecta, whereas the scattering in other lines 
(e.g., Fc ll) is neglected. The quasi-continuum intensity, the 
slope of the emergent spectrum, and the relative line inten- 
sity are free parameters. 

The calculated Ca ii doublet assuming a constant ve- 
locity and an occupation optical depth Qo = (model DCO) 
and Qo — 1 (model DCl), and assuming a linear velocity dis- 
tribution with the same Qo values (models DLO and DLl, 
respectively) are shown in Fig. 7 overplotted on the spec- 
trum on day 244 (Wang et al. 2004) . The computed spectra 
are smoothed by a gaussian filter with FWHM=7 A to take 
into account the finite spectral resolution. Velocities Vs and 
Van of the standard interaction model as well as other pa- 
rameters are presented in Table 1. The table gives also Qo, 
the optimal cut-off radius (rc) and the optical depth of the 
SN ejecta in the Ca ll 3934 A fine (rsn). We found that ef- 
ficient scattering in the SN ejecta is needed with Tsn > 3 in 
Ca II 3934 A. The optimal cut-off radius for all the models, 
rc = 0.95, coincides with the inner boundary of the line- 
emitting shell. 

All four models shown in Fig. 7 reproduce the major 
properties of the observed Ca ll doublet, particularly, the 
width, shape and the skewing towards the blue. The de- 
viation in the red wing may be related to the behavior of 
the underlying quasi-continuum in the observed spectrum. 
Models with constant and linear velocity laws in the shell are 
practically indistinguishable. We will use below only the lin- 
ear law model that presumably mimics the possible random 
component of the velocity field in the shell. 

Using models DLO and DLl, we also simulated the Ca ii 
triplet (Fig. 8). The model profiles are overplotted on the ob- 
served profile in the Subaru spectrum on day 234 (Deng et 
al. 2004) . Parameters of the corresponding models (TLO and 
TLl) including velocities of the standard model are given in 
Table 1 together with the optimal contribution of the O I 
8446 A line to the 8500 A feature C(0 l). The overall fit of 
the model TLO is satisfactory. The derived contribution of 
the O I 8446 A fine to the 8500 A feature C(0 i) = 0.25 is 
constrained with an uncertainty of ±0.05. However, model 
TLl with Qo = 1 shows a pronounced central absorption 
feature caused by the resonance non-local absorption in the 
shell. This feature is inconsistent with the observed profile. 
We conclude that the fragmented shell should be character- 
ized by a low probability of non-local radiation coupling in 
the Ca II triplet. We will return to this issue below (Sec. 
3.3). 

The observed 8500 A profile is slightly broader than the 
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Figure 7. Calculated profiles of the Ca II 3934, 3968 A doublet 
(thick line) overplotted on the spectrum on day 244 (Wang et 
al. 2004). Upper panels show models with constant velocity in 
the shell and different Qo (model DCO, left, and model DCl, 
right). The lower panels are the same except for the linear velocity 
distribution in the shell (model DLO, left, and model DLl, right). 
Two vertical bars show the rest frame positions of 3934 A and 
3968 A lines. 

model feature, especially in the wings. However, one rimst 
be cautious about attributing the wings to the Ca ll triplet 
because the quasi-continuum also contributes to the 8500 
A feature, as will be demonstrated below. Given the un- 
certainties of the behavior of the quasi-continuum and the 
satisfactory fit of the Ca ll doublet and triplet, we conclude 
that the Ca II features are consistent with a spherical model 
and expansion velocities predicted by the interaction model. 

We performed computations with velocities multiplied 
by a factor larger than unity. We found that slightly larger 
expansion velocities (by ^ 5%) produce better fits to the 
8500 A features. However, this conclusion is uncertain be- 
cause of a possible contribution from the quasi-continuum. 

3.3 CDS temperature and S/So ratio 

The contribution of the O I 8446 A line to the 8500 A fea- 
ture found on day 234 (~ 25%), combined with the total 
flux « 4.7 X 10"^* erg s"^ cm"^ (Deng et al. 2004), im- 
plies a flux in the O I 8446 A fine « 1.2 x 10"" erg s"^ 
cm~'^, and « 3.5 x 10"^'' erg s~^ cm~^ in Ca ll triplet. 
Given the Ca ll doublet flux ~ 10^^'' erg s^^ cm^^ in the 
same spectrum (Deng et al 2004), one gets the triplet-to- 
doublet ratio F(8579)/F(3945) « 3.5. Assuming that both 
Ca II excited terms (^D and ^P) have the same excitation 
temperature and using the scaling F oc uB^{T) from the 
ratio F(8579)/F(3945) « 3.5 we obtain T « 4430 K. 

The absolute luminosities of the Ca ii features can be 
used to estimate T and the area ratio S/So simultaneously 
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Figure 8. Calculated profile of the Ca II infrared triplet blended 
with O I 8446 A for the models TLO (left) and TLl (right) over- 
plotted on the observed profile on day 234 (Deng et al. 2004). 



(see equation [4]). To make the procedure reliable, we aban- 
doned the assumption that the excitation temperature for 
both excited terms is the same. Instead, the statistical equi- 
librium for a three level Ca ll model is solved with coUisional 
and radiative transitions taken into account. The radiation 
transfer in a fragment is treated in the static average escape 
probability approximation. The adopted Ca abundance is 
0.03 by mass, the same as used in Section 3.1.1. The ion- 
ization fraction of Ca ll is calculated in the LTE approxi- 
mation, which implies that Ca ii is the dominant stage for 
a wide range of temperature, 4000 — 6000 K. The solution 
of statistical equilibrium then yields the excitation temper- 
ature in both Ca ll transitions as a function of the electron 
temperature Tc and S/So, which permits us to calculate the 
luminosities of the Ca ii features (equation [4]) and to de- 
termine Te and S/So from a comparison with the observed 
luminosities. 

The Ca ii levels are found to be well thermalized, so 
to reproduce the observed luminosities of Ca ii features we 
need only two parameters: Te and S/So, assuming ^ = 0.05. 
Figure 9 shows the behavior of the luminosities of the Ca ii 
features as a function of S/So for the optimal temperature 
Te = 4430 K. The plot clearly demonstrates that the ob- 
served luminosities of the Ca ii features can only be repro- 
duced for a large area ratio, S/So ~ 10^. This result pro- 
vides strong evidence that the dense material of the shocked 
ejecta experienced substantial deformation, fragmentation 
and mixing in the intershock layer. 

We studied the sensitivity of the luminosity of Ca ii lines 
to the adopted mass locked in line-emitting fragments (Mc) 
and found that in the range of S/ So < 300 the luminosity of 
both Ca II multiplets are well reproduced for Mo > 10~^ Mq 
and cannot be reproduced for smaller mass, Me < 10"^ Mq. 
This indicates that the amount of the Ca-rich matter locked 
in the Ca ll line-emitting gas is > 0.01 Mq (for the mass 
fraction of Ca ^ 0.03). 

The conclusion that the area ratio is as large as ~ 10^ 
raises a serious problem. The local occultation optical depth 
of the shell with an area ratio S/So ~ 10^ is Qo ~ 1- This 
seems inconsistent with the claim (Sec. 3.2) that the non- 
local radiative coupling in the Ca ii triplet is weak. For the 
spherical shell, this requirement implies Qo <S 1- In an at- 
tempt to resolve the inconsistency, we considered a model of 
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Figure 9. Calculated luminosities of the Ca ii doublet and triplet 
on day 234 as a function of the area ratio for Te = 4430 K with 
zero albedo {thick line) and finite albedo (dashed line) of dense 
fragments. The observed luminosities of the Ca li doublet and 
infrared triplet are plotted by dotted lines. Subscripts T and D 
refer to triplet and doublet feature. 



the Ca II feature formation with a finite albedo lj of dense 
fragments for the incident line photons. To determine lu we 
performed a Monte-Carlo simulation of the history of the 
incident Ca ii doublet and triplet photons assuming a three 
level model of Ca ll. The line luminosities for the finite 
albedo are shown in Figure 9. The required value of the 
area ratio in the new model is only slightly lower compared 
with the previous model, so the problem of large Qo remains 
unsettled. 

A way out might be in the assumption that the accepted 
model of the fragmented shell is too simple to predict the 
upper limit for the occultation optical depth from the Ca ii 
triplet. For instance, the large scale corrugation of the shell 
and/or patchy structure of the line-emitting zones might re- 
sult in voids lacking the absorbing dense gas at the resonance 
points. As a result, the large local value of Qo may be recon- 
ciled with the low absorption probability for the non-local 
interaction in the Ca ii triplet. 

3.4 Quasi-continuum model 

We now check our conjecture that the quasi-continuum of 
SN 2002ic is produced by the emission of the fragmented 
cool shocked ejecta in numerous lines of Fe-peaJj elements. 
To this end, we compute the emergent spectrum assuming 
a model of the fragmented shell width AR/R = 0.1 for two 
values of the area ratio S/So = 50 and S/So = 1- The shell 
patchiness indicated by the Ca ll triplet modelling will be 
described by the surface (and volume) filling factor /qc of 
quasi-continuum emitters. It should be emphasized that the 
filling factor refers only to the volume fraction of mixture 
and not to the volume fraction of the dense material that 
is, of course, much lower. The radiation transfer in the shell 
will be computed assuming a spherical shell, with two cor- 
rections. First, the emitting volume will be multiplied by the 
factor /qc, and, second, we take into account a finite proba- 
bility of the photon escape through the holes of the opposite 
side of shell, w = |/i|(l — /qc)- Here /i, the cosine between 
the photon wave vector and the radius, describes approxi- 
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mately the angular dependence of w. The unshocked ejecta 
wiU be considered as transparent and the net emission of 
this ejecta will be neglected. We also omit the presence of a 
shocked H-rich cool gas in the intershock layer. 

We consider the seven most abundant Fc-pcak elements, 
Ti, V, Cr, Mn, Fe, Co, Ni, assuming their relative abun- 
dance to be solar. The latter is a sensible approximation for 
the explosive nucleosynthesis in the full incineration regime. 
Doublet and triplet lines of Ca ll are also included in the 
line list. The initial line list is extracted from the VALD 
data base (Kupka et al. 1999) assuming conditions in a red 
giant atmosphere with T^g = 5000 K and gravitational ac- 
celeration log g = 1.5. The list contains about 19000 lines of 
neutral and singly ionized species with a line optical depth 
> 10"* in the range of 3500 - 10000 A. Using this list, we 
compile the input line list for a particular density, assum- 
ing Saha ionization with a temperature Tion and Boltzmann 
excitation with a temperature Tox. The line optical depth is 
restricted by the condition r > 10^*. For the relevant condi- 
tions, the most abundant species are singly ionized (02 ~ 1), 
while the ionization fraction of neutrals is 0i ~ 10~^ and 
that of double ions is 03 ~ 10~®. 

The radiation transfer is calculated using a Monte Carlo 
technique. The photon interaction with a dense fragment is 
treated as a partial absorption with an albedo w. We adopt 
u! = 0.1 for all the lines, although we found that the re- 
sult is practically the same in the case of pure absorption 
{uj — 0). Scattering is assumed to be isotropic for the back 
and forward scattering with a diffusion reflection probabil- 
ity r/(l + t). The computed spectrum is reddened using 
E{B -V)= 0.06 (Deng ct al. 2004). 

An extensive search in the parameter space (Tion, ?ox, 
/qc) for the model with S/So = 50 led us to conclude that a 
best fit of the quasi-continuum in the spectrum on day 244 
is attained for the parameters Tion = 6000 K, Tex = 4300 K, 
and /qc = 0.125 (Fig. 10a). Note, the result is not very much 
sensitive to the area ratio; we adopt S/So = 50, not 100 that 
follows from Ca n line analysis, because the latter value leads 
to a too large occultation optical depth (see discussion in 
Section 3.3). On the other hand, the optimal model with 
the small area ratio, S/So = 1 and /qc = 1 (Fig. 10b) seems 
to be inconsistent with the observed shape of the quasi- 
continuum. This implies that the model with a large area 
ratio is preferred thus supporting the conclusion made from 
Ca II line analysis. 

The fit is good despite the fact that the quasi-continuum 
model is rudimentary; the number of free parameters is 
small. The agreement suggests that we are on the right 
track in the interpretation of the quasi-continuum. This re- 
sult makes the identification of [O l], [O ll], [Ca ll] lines for 
the 6400 A and 7300 A bumps (cf. Deng et al. 2004) un- 
neccessary. Otherwise, there would be serious problems in 
the interpretation of the identified emission lines. Another 
remarkable fact is that the quasi-continuum contributes to 
the 8500 A emission feature thus alleviating the problem of 
the origin of broad wings in the Ca ii/ O I 8500 A feature. 
We checked the contribution of different ions to the quasi- 
continuum and found that, as expected, the Fe ii lines are 
dominant. This elucidates the role of the ionization temper- 
ature as a fitting parameter that regulates the contribution 
of neutral species, Fe I in particular. The quasi-continuum 
in the region A < 5500 A is very reminiscent to that in early 
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Figure 10. Synthetic spectra for the model of a fragmented shell 

ovciplottcd on the observed spectrum of SN 2002ic on day 244 
(Wang et al. 2004). Panel a shows the model with a large area 
ratio {S/So = 50, see text), whereas panel b shows a model with 
S/So = 1. Panel c shows the same model as in panel a but with 
the addition of a Ca o line-emitting zone. 



SN 1988Z, which also has been identified with Fe ii lines 

(Clmgai & Danzigcr 1994). 

Interestingly, the smooth quasi-continuum in the range 
A < 4500 A (Fig. 10a) is the result of radiation saturation at 
the Planck intensity due to the high density of lines in this 
band that results in the large optical depth of the shell, viz. 
Tshcii ^ {Q)Nunc > 1, whore Nu^c is the number of linos in 
the characteristic wavelength interval \{vs/c){AR/R) with 
the average local occultation optical depth (Q). The sat- 
uration efi^ect is absent in a model with small area ratio 
S/So ~ 1 (Fig. 10b). Because the line density in the ul- 
traviolet (UV) region is very high, wc expect that the IJV 
quasi-continuum should retain a smooth black-body appear- 
ance with T ~ 4300 K in the short wavelength region 
(A < 3500 A). 

As shown in Fig. 10a, the quasi-continuum model is 
unable to produce strong Ca II emission features. This is re- 
lated to the fact that the quasi-continuum is a superposition 
of a large number of weak lines afi^ected additionally by the 
line overlapping. We propose, therefore, that the Ca ll emis- 
sion and quasi-continuum should form in different zones of 
the shocked ejecta. These zones, presumably, correspond to 
Fe-poor (Ca-rich) gas and Fe-rich gas, respectively. The ob- 
vious physical reason for these distinctive zones is the abun- 
dance stratification of SNc la with the inner layers incin- 
erated up to Fe-peak elements and the outer layers having 
undergone incomplete burning. In this situation one expects 
that the outer Fe-poor gas that has passed through the re- 
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verse shock first has enough time to spread over tiic large 
volume of the intershock layer, whereas the inner Fe-rich 
material shocked recently occupies only a small fraction of 
the intershock layer. 

To illustrate the proposed scenario for the Ca ll lines, 
we consider a combined model of the fragmented shell: apart 
from the quasi-continuum line-emitting zones suggested by 
the model with /qc = 0.125, we consider the Ca II line- 
emitting zones with a filling factor C(l — /qc), where C is 
a fitting parameter. Wc assume a Ca abundance ~ 0.03 by 
mass (the result only weakly sensitive to Ca abundance) and 
neglect the Fe-peak elements in the Ca ii line-emitting zone. 
For Ca ll lines, the scattering in the unshockcd SN cjocta 
on Ca II is taken into account as described in the previous 
section. The same area ratio S/So — 50 and Tex = 4300 
K are adopted in the Ca ii and quasi-continuum zones. We 
omit here the contribution of O I 8446 A. 

The result presented in Fig. 10c is obtained for C = 0.8, 
which suggests that ~ 70% of the fragmented shell volume 
is occupied by Fe-poor Ca ll line-emitting material. The fill- 
ing factor again refers only to the volume fraction of the 
mixture and not the dense material. Although we do not 
intend to describe the Ca ll profiles and intensities in detail, 
the combined model provides a satisfactory representation 
of the overall spectrum, except for Ha and H/3. The mod- 
elling thus confirms the proposed scenario based upon the 
expected abundance stratification in the outer layers of the 
SN ejecta. 

We now address the question of how low the abun- 
dance of Fe-peak elements in the Ca ii line-emitting mat- 
ter can be and retain sufficient contrast between the quasi- 
continuum and the Ca II lines, especially the Ca II dou- 
blet that resides in the line-rich spectral range. The an- 
swer depends on the assumed mass locked in the Ca ii line- 
emitting dense gas (Me). We adopt Me m 0.02 Mq and a 
Ca abundance 0.03, which is close to the admissible lower 
limit (Mo > 0.02 Mq) and is consistent with the mass and 
composition of the Fe-poor Ca-rich layer produced by the 
delayed detonation model (Iwamoto et al. 1999). For these 
parameters and the values of Tion and Tex adopted above 
for the case of S/So = 50, we found that the mass fraction 
of Fe-peak elements should not exceed 2 x 10~* or ~ 0.1 of 
solar abundance. This upper limit indicates a low primordial 
metallicity of the pre-SN and, possibly, of the host galaxy 
as well. 

To summarize, the major observational properties of 
the SN 2002ic spectrum at an age of ~ 240 days may be 
explained in the framework of a model of a spherical, locally 
inhomogeneous shell with an expansion velocity predicted 
by the standard hydrodyriamic model. 



4 DISCUSSION AND CONCLUSIONS 

The aim of this paper was to develop a unified model for 
the light curve and spectrum of SN 2002ic assuming ap- 
proximate spherical symmetry. Our principal result is that 
the Ca II line profiles and the quasi-continuum are consis- 
tent with the interaction model provided that the SN la has 
the maximum possible kinetic energy. 

The high explosion energy is in line with the high peak 
absolute magnitude of SN 2002ic. According to Hamuy et 



al. (2003), at light maximum SN 2002ic was as bright as 
M(V) f» -20.3 {Ho = 65 km s"^ Mpc"^). Given the con- 
tribution of the veiling continuum of about 0.3 mag on Dec. 
SinV band (5000 - 6000 A) (Hamuy et al. 2003), we con- 
clude that the peak absolute magnitude of the SN la was 
M{V) « -20, close to that of extremely bright SN la, 
SN 1991T (M(y) w -20.2, cf. Fisher et al. 1999). The mag- 
nitude at the maximum coincides also with that predicted 
by the detonation model DETl (Hoflich & Khokhlov 1996). 
This coincidence is another justification for our choice of 
the high kinetic energy for the SN 2002ic event. Yet, as we 
noted above, the strong Ca ii emission lines indicate that 
the explosion should be different from the DETl model in 
the sense that the degree of incineration in the outer layers 
should be lower to produce enough (> 10~^ Mq) Ca-rich 
Fe-poor material. 

Modeling of the spectra of the fragmented cool dense 
shell formed in the intershock layer provides strong evidence 
that the quasi-continuum and Ca ll lines are emitted by 
different zones. This presumably refiects the stratification 
of the composition of the SN ejecta in the sense that the 
outer layers responsible for the Ca ii lines are Fe-poor, while 
the inner layers responsible for the quasi-continuum are Fe- 
rich. The Fe-poor ejecta are shocked earlier and, therefore, 
spread throughout the intershock layer, while Fe-rich ejecta 
are shocked later and thus have not yet mixed with the Fe- 
poor material at t ^ 240 d. This picture predicts that at a 
later epoch (t > 240 d) both Fe-poor and Fe-rich shocked 
components should be mixed with each other. As a result, 
the Ca II doublet must disappear at a later epoch due to the 
mixing of the Ca ll line-emitting gas with optically thick 
quasi-continuum emitters. SN 1999E showed a gradual dis- 
appearance of the Ca ll 3945 A doublet in spectra between 
days 139 and 361 (Rigon et al. 2003), in accord with these 
expectations. 

The high luminosity of the Ca ll doublet and triplet 
in SN 2002ic on day 234 indicates that the dense shocked 
ejecta form a structure with a large area ratio, S/So ~ 10^. 
A high area ratio {S/So ~ 50) is also preferred by the quasi- 
continuum model. The area ratio, which is common in the 
treatment of the physics of mixing fiuids (e.g., Catrakis et 
al. 2002), hcis not been used before in the interpretation 
of optically thick emission lines from an intershock layer. 
It was recognized earlier that the interpretation of the X- 
ray and radio emission of an interacting SN depends on the 
presence of a mixing layer, while the Ha emitted by the 
mixed dense material may reveal the fine structure related 
to the RT spikes (Chevalier & Blondin 1995). Here we add 
to the list of effects of the mixing layer the crucial role of 
the large area ratio for the interpretation of the Ca ll lines 
and quasi-continuum in SN 2002ic. 

The fact that the area ratio in SN 2002ic on day 234 
is large, S/So ~ 10^, suggests that we observe an advanced 
stage of the mixing of dense shocked ejecta in the inter- 
shock zone. Experiments on the mixing of fluids show that 
the interface between fluids reveals fractal behavior. Specif- 
ically, the cumulative area of the interface above the scale 
A is 5 oc A~™, where a « 0.35 (Sreenivasan et al. 1989). If 
this law also holds in the case of the mixing of the dense 
shocked ejecta in the intershock layer of SN 2002ic, then the 
dense gas should produce structures with scales as small as 
X/R DC {S/So)'^^" ~ 10"^. The minimum Unear scale that 



© 0000 RAS, MNRAS 000, 000-000 



Circumstellar interaction of SN 2002ic 11 



corresponds to the thickness of the dense fragments should 
be even smaher, ~ 5{So/S) ^ 10~^R, where 5 is the origi- 
nal thickness of the CDS (Section 3.1). It would be interest- 
ing to check the possibility of creating the mixing layer in 
SN 2002ic with S/So ^ 10^ in numerical experiments using 
high resolution 3D hydrodynamical simulations, although 
it seems to be beyond the capabilities of present-day com- 
puters. To estimate the area ratio from a numerical experi- 
ment, one needs to find the average number of intersections 
of dense peaks with the radius in the intershock layer TVis- 
The area ratio is then S'/S'o « 2A^is. 

We do not yet fully understand all the mechanisms re- 
sponsible for the major spectral features seen in the late 
emission from SN 2002ic. In addition to the Ha lino, an im- 
portant issue ignored in our study is the problem of the O I 
8446 A emission, although it is not a dominant constituent 
of the 8500 A bump anymore. The oxygen excitation mech- 
anism should predict the absence of strong O I 7774 A in 
SN 2002ic. The fluorescence pumping of O l 8446 A in the 
Ha-emitting gas (Deng et al. 2004) meets this requirement. 
This mechanism predicts a comparable intensity of O 1 11287 
A. The absence of this line in the spectrum of SN 1999E 
around ~ 200 d (Rigon ct al. 2003, their Fig. 6) casts doubt 
on the fluorescence mechanism. An alternative possibility 
is that the O I 8446 A line is emitted by the oxygen-rich 
matter in the outermost layers of the SN ejecta due to non- 
thermal collisional excitation by secondary electrons that 
are created by high energy photoelectrons. Excitation by 
secondary electrons also yields a high O I 8446 A/O I 7774 
A intensity ratio (Stolarski & Green 1967). An interesting 
possibility that Lya pumping line of Fc ll 8451 A might con- 
tribute to 8500 A band instead of O I 8446 A, as was the 
case, possibly, for SN 1995N (Fransson et al. 2000), seems 
unlikely here because of the absence of another strong Lya 
fluorescent band of Fe ii at ~ 9150 A. Despite of a possible 
presence of a large amount of oxygen in the cool shocked SN 
gas (up to ~ 0.1 Mq), the strong [O l] line is not expected. 
Adopting Boltzmami population for T — 4400 K and oxygen 
mass of 0.1 Mq one gets the luminosity of this line of about 
^ 10^^ erg s~^, i.e., a factor of ~ 10~^ lower compared to 
the luminosity of the observed "6300 A" emission band on 
day 244. 

In view of the possibility that our hydrodynamic model 
may underestimate the expansion velocity of the SN/CSM 
interface by about 5%, we discuss the implications for our 
model. The emergent X-ray luminosity from the forward 
shock in the standard model is comparable or even exceeds 
the optical bolometric luminosity after maximum light (see 
Fig. 1 and Fig. 3). Hence, there is an extra energy reservoir 
which might, in principal, be transformed into the optical 
band under the proper conditions. This possibility could be 
realized in a model of interaction with a clumpy CSM. As 
a result of (i) softer X-rays from cloud shocks than for a 
smooth CSM, and (ii) the larger geometrical probability of 
absorption of X-rays from clouds penetrating the SN ejecta, 
the optical output from the interaction should be larger. 
Therefore, the same bolometric luminosity can be produced 
for a lower average CS density and higher expansion velocity 
of the SN/CSM interface compared with a smooth CSM. 

The spherical geometry of the CS interaction of 
SN 2002ic, advocated here, finds additional support from 
statistical arguments already mentioned in the Introduction. 



We know three events of SNe la in a dense CSM (SN 1997cy, 
SN 1999E, SN 2002ic) and they show similar widths and 
shapes of spectral features at a similar phase (see Deng et 
al. 2004, their Fig. 1). This fact suggests that a bi-polar 
structure of the SN/CSM interface with high aspect ratio is 
unlikely. 

At first glance the picture of spherical interaction is 
incompatible with the detection of intrinsic polarization of 
SN 2002ic radiation at a level of - 0.8 % (Wang et al. 2004). 
Simple estimates show that the Thomson scattering in a 
moderately aspherical CS envelope with tt 0.1 is not 
able to produce such a high polarization. The issue may be 
resolved in a model in which the polarization is acquired in 
the scattering of the SN radiation on the dust in a distant 
(> lO^'' cm) aspherical CSM. 

The resemblance of the light curves for all three 
SN 2002ic-like events (see Deng et al. 2004) suggests a simi- 
lar regime of mass-loss from their prcsupcrnovac and a com- 
mon origin of the progenitors. As mentioned before (Chugai 
& Yungelson 2004), the uniformity lends support for the 
SN1.5 evolutionary scenario suggested by Iben & Renzini 
(1983). Yet we do not rule out the binary scenario. In this 
case, the WD companion should be a supergiant with the 
mciss > 2 Mq (the CS envelope, « 1.6 Mq, plus the super- 
giant core, > 0.5 Mq). In the binary scenario, there is no 
clear reason for the supergiant to have experienced strong 
mass loss just before the SN explosion. The synchronization 
of the vigorous mass loss and the SN event might be a more 
natural outcome in the single AGE star scenario with the 
C/O core approaching the Chandrasekhar limit. 

A general belief is that Chandrasekhar mass white 
dwarfs cannot be produced from intermediate-mass single 
stars because the growth of the C/O core during the thermal 
pulsation stage is prevented by the removal of the hydrogen 
envelope (Blocker 1995, and references therein). In this re- 
gard, it is highly remarkable that at least two supernovae 
(SN 2002ic and SN 1997cy) exploded in dwarf galaxies, as 
already noted by Chugai & Yungelson (2004). Dwarf galax- 
ies are characterized by low metallicity, which thus might 
explain why the CO core of an asymptotic giant branch 
star is successful in attaining the Chandrasekhar limit before 
the hydrogen envelope is lost (Chugai & Yungelson 2004). 
An evidence for the low metallicity of the presupernova of 
SN 2002ic comes from our modelling of the Ca II lines and 
quasi-continuum (Sec. 3.4). The possibility of SN1.5 explo- 
sions in metal-poor galaxies was suggested also by Zijlstra 
(2004). It would be interesting to check whether all the host 
galaxies of SN 2002ic-like supernovae are actually metal de- 
ficient. If the suggested link between SN 2002ic-likc events 
and low metallicity is correct, then their occurrence rate 
relative to normal SN la must increase with the redshift for 
z > 2 — 3, when the average cosmic metallicity was lower 
than solar. 

If SN 2002ic-likc events stem from SN 1.5 explosions, 
their low fraction of SN la, < 1%, can be used to estimate 
the the fraction of white dwarfs created with a mass sim- 
ilar to the Chandrasekhar mass. If only ~ 1% of all WDs 
can become SN la, this suggests that only < 10~* of white 
dwarfs in our Galaxy were formed with a mass close to the 
Chandrasekhar mass. Since these white dwarfs are likely to 
have formed in a metal poor environment, they are prob- 
ably very old (age ~ 10^° yr) and thus of low luminosity. 
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< 3 X lO"'^ Lq (Prada Moroni & Straniero 2003). The num- 
ber of known cool white dwarfs is low, < 40 (Salim et al. 
2004), and will probably remain too low for the nearest fu- 
ture to test the idea whether high mass white dwarfs were 
produced during early epochs when the mctallicity was low. 

The possible link between SN 2002ic-like supernovae 
and hypothetical SN 1.5 events, combined with arguments 
in favour of a high energy SN la event in this case, raises 
the question of whether conditions in a hot Chandrasekhar 
mass C/O WD that just got rid of the remnants of the hy- 
drogen envelope favour a WD explosion with a high kinetic 
energy of eject a. 

The SN 1.5 origin of SN 2002ic suggests that the initial 
mass at the main sequence is ~ 8 Mq. With the Chan- 
drasekhar core and 1.6 Mq in the dense CS envelope we 
expect that about 5 Mq had to be lost at the previous 
AGB stage. This means that outside the dense CS envelope 
(r > 3 X 10^^ cm) there should be a dense red supergiant 
wind with a velocity of the order ~ 10 km s~^ and a mass 
loss rate of ~ 10~'' Mq yr~^. Keeping in mind the analogy 
with SN 1979C, this wind may be revealed in the infrared 
(IR) as an IR echo (Dwek 1983), and in the radio as a result 
of CS interaction (Chevalier 1982a, Woilor ot al. 1986). How- 
ever, the radio emission from the distant SN 2002ic 285 
Mpc) may pose a difficult challenge to detect. If the radio 
luminosity is comparable to that of SN 1979C (Weiler et al. 
1986) then the expected radio flux from SN 2002ic around 
1000 day should be of the order of ~ 0.03 mjy in the 6 — 20 
cm band, which is close to the detection limit of VLA. A 
more promising target for the detection of the late time in- 
teraction of SN 2002ic-likc supernovae with the presumed 
AGB wind is SN 1999E at a distance of ~ 112 Mpc, or any 
closer SN 2002ic-like event that will be discovered in the 
future. Note, that at the epoch of the interaction with the 
dense CS shell {t < 500 d), the radio emission is severely 
attenuated by free-free absorption at wavelengths A > 1 
cm. With the parameters used in our models, and adopting 
temperatures for the shell and wind in accordance with the 
models of Lundqvist & Fransson (1988), we find that free- 
free absorption is a likely reeison why the observations of 
Bcrger & Sodcrbcrg (2003) and Stockdale et al. (2003) did 
not detect radio emission from SN 2002ic. 

A better understanding of the phenomenon of 
SN 2002ic-like events requires an accurate energy audit, for 
which the direct detection of the X-ray emission as well as 
IR measurements arc important. SN 2002ic-likc events at 
an age of 200 — 500 days are excellent targets for Chandra, 
XMM and INTEGRAL. 
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